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Abstract. Using exact relations between velocity structure functions [TJ [3] and 
neglecting pressure contributions in a first approximation, we obtain a closed system 
and derive simple order-dependent rescaling relationships between longitudinal and 
transverse structure functions. By means of numerical data with turbulent Reynolds 
numbers ranging from = 320 to 5Ra = 730, we establish a clear correspondence 
between their respective scaling range, while confirming that their scaling exponents 
do differ. This difference does not seem to depend on Reynolds number. Making use of 
the Mellin transform, we further map longitudinal to (rescaled) transverse probability 
density functions. 
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1. Introduction 

Intermittency is an ubiquitous feature of fluid turbulence: the scaling properties of 
flow quantities differ from Kolmogorov's mean field theory [3J. For instance, in the 
inertial range of scales where flow properties are assumed to be independent of the 
details of energy injection and dissipation, the velocity increments do not have a 
monofractal structure. In homogeneous isotropic turbulence (HIT), two directions only 
matter in the computation of a velocity increment: the longitudinal one, A r u taken 
along the separation and the transverse one, A r v in which the difference of velocities 
components perpendicular to the separation are computed. Velocity structure functions 
are then defined as the ensemble average: S n>m (r) = ((A r u) n (A r v) m ). For HIT, the 
von Karman-Howarth relationship taken in the inertia range, leads to Kolmogorov's 
4/5th law: (Aw r ) 3 = — |(e)r, where (e) is the mean energy dissipation rate per unit 
mass. While a monofractal inertial range behavior would then lead to S n> o oc r^ n 
with ( n ~ rt/3, intermittency means that ( n is a non-linear (concave) function of n. 
Numerous works have been devoted to the study of the functional form of ( n . We 
focus on the possible link between the longitudinal and transverse structure functions 
S n ,o(r) and So iTl (r). There exist theoretical arguments that longitudinal and transverse 
show the same scaling [5]. However, both experimental data [6J IH1 [9] and numerical 
simulations pHl EH EGJ [13] show consistently different scaling exponents for longitudinal 
and transverse structure functions. Whether this difference can be attributed to a 
persistent small scale anisotropy [TJJ [15] or to a finite Reynolds number effect p] is 
an unsolved question to which we will come back below. Here, we note that in the 
case of the direct cascade in electron-magnetohydrodynamic [16] it was demonstrated 
numerically that the difference vanishes with increasing numerical resolution and thus 
is a finite Reynolds number effect. 

In this article, we rather focus on the correspondence between scaling ranges of 
the longitudinal and transverse structure functions. Our approach is based on the 
observation that even though the (real space) velocity field of a turbulent flow coarse- 
grained at a scale r is not smooth, the structure functions are smooth ( different iable) 
functions of r. We thus use the structure of the Navier-Stokes equation together 
with assumptions to derive constitutive relationships between S n ,o(r) and So, n (?")- 
Specifically, we shall neglect the contributions from the pressure term. We start with 
exact scaling expressions derived by Hill [T], Hill and Boratav [2] and Yakhot [3]; 
we then obtain rescaling relationships between longitudinal and transverse structure 
functions. Our first finding is that, after rescaling, the longitudinal and transverse 
structure functions share the same inertial range, i.e. the same width in the extent of 
of scales where self-similarity is observed. This is important because the question of 
the location and span of the inertial range is often an issue in the analysis of turbulent 
data at (necessarily) finite Reynolds number. We stress, however, that the value of 
the longitudinal and transverse scaling exponents do differ. A second outcome of our 
simple ansatz is a direct mapping, using the Mellin transform, of the transverse and 
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longitudinal probability density functions (PDFs). Differences which persist after the 
mapping are then due to the effect of the neglected terms, as pointed in some previous 
attempts by Yakhot [3j and Gotoh and Nakano p2] . 



2. Rescaling relations between longitudinal and transverse structure 
functions 



Our calculation trace back to the observation by Siefert and Peinke [18] that the 
Karman equation (see Karman and Howarth [19]) relating second order longitudinal 
and transverse structure functions can be interpreted as a Taylor expansion of a smooth 
function. To see this, we start with the Karman equation 

SoM = S 2 ,o(r) + —S^r) , (1) 

which is exact, and contains no contribution from the pressure - it is a statement 
of incompressibility. Siefert and Peinke [18] observed that the structure function is a 
smooth function of r and that if the scale r is chosen in the inertial range i.e. "small" 
compared to the integral scale L, eqn. (P) can be seen as a Taylor expansion: 

SoAr) « S 2i0 (r + "-) = S 2fi (^r) , (2) 

where the function 5*2,0 is expanded about r for consistency with the exact 
relationship flTJ). In [18J evidence from experimental data Taylor-based Reynolds 
numbers between 180 and 550 was presented to support this view. The success 
of the approach introduced by Siefert and Peinke [H] motivated us to extend their 
reinterpretation of differential relations to structure functions of higher orders, making 
use of the exact relationships derived by Hill & Boratav [2], Hill pQ and Yakhot [3]. Hill 
derived these relations directly by inventing a clever matrix algorithm which allowed 
him to efficiently simplify the derivation and calculations. Yakhot [3], on the other 
hand, derived an equation for the characteristic function Z — (A • Au r ) where Au r 
denotes a velocity increment over the distance r. Structure function relations can then 
be obtained by differentiating the characteristic function Z. 

As an illustrative example, consider the relation for even order mixed structure 
functions derived by Yakhot [3] : 

dS 2nfl ^ 2 2(2n -1) _ 

7; 1 >->2n,0 H J2n-2,2 — L/„ + G / . [6 

The term C p contains contributions from the (unknown) pressure field and is the reason 
why the system cannot be closed. The term Cf contains contributions from the large 
scale forcing and can safely be ignored in the inertial range as proposed by Kurien 
and Sreenivasan [20]. These authors also analyzed and compared these relations to 
measurements in atmospheric turbulence at a Taylor-based Reynolds-number of about 
10 700. One of their findings was that for even order structure functions the pressure 
contributions can be an order of magnitude smaller than the terms directly related to 
the structure functions. A detailed numerical study on the role of the pressure term 
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has been realized by Gotoh and Nakano [TTJ . In order to obtain closed expressions, we 
shall hereafter neglect the pressure contributions. Although this assumption is quite 
crude, we are already able to capture the dominant features of the relationship between 
longitudinal and transverse structure functions, such as amplitudes and common inertial 
range. Pressure (or energy injection) contributions will then appear as departures from 
predictions of this closed system of equations. 

In order to demonstrate the procedure, we start with formulas for the 4-th order 
structure functions and neglect contributions from the pressure and the large scale 
forcing (see also eqn.(ll) and (13) in [20J) 

3S 2i2 (r) « S 4)0 (r) + —S^r) « S 4 , (^r) , (4) 

^S 0A (r) « S 2>2 (r) + r -^S 2 ,2{r) « S 2 , 2 (^r) , (5) 
which can be combined into 

^4,0 (||r) « So A (r) ■ (6) 

For the 6th order structure functions we get similarly (see eqn.(12), (15) and (14) in |20j) 

v d 3 
5S 4 , 2 (r) « S 6>0 (r) + -^-S 6fi (r) « S 6>0 (-r) , 
2 or 2 

S 2A (r) « Sw(r) + —^(r) « S 4 , 2 (^r) , 

^ 0i6 (r) « 5 2)4 (r) + ~S 2 4r) « 5 2 , 4 (^r) . 
5 6 ar 6 

Again, combining these equations results in the simple relation 

^6,0 Q^ r J ~ S oA r ) ■ ( 7 ) 
In general, the rescaling for even order structure functions reads 

Sn,o • • • ^r) = S n , ( 2n ^ 2 2) +1) r) - S 0>n (r) (8) 

In order to demonstrate that the Taylor expansion is valid also for higher order 
structure functions, we look at the differential relation for the 4th order structure 
function obtained from eqn. (j4]) and (J5]) 

c ( r] _a ( r \ | 7 ^4o(r) l ^ 2 »S 4 o(r) 2 

5 4(r)-5 40 (r) + -^— r + -^-^r . (9) 

In Fig. [T]we compare the longitudinal (black), transverse (blue), rescaled longitudinal 
(red) and the one using the differential relation (green). The difference between the 
rescaled longitudinal and the one using the differential relation is negligible. 

In order to test our results, we use numerical data from pseudo-spectral simulations 
of incompressible Navier-Stokes turbulence, as described in |2T] - the LaTu code. A 
statistically stationary flow is maintained by keeping constant the Fourier modes in the 
two lowest shells. All results are averaged over several large-eddy turn-over times (over 
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«rm S £k V_ dx 7/ Tj, L T L N 3 

730 0.192 3.8 • 10~ 3 1. • KT 5 1.53 • 10~ 3 7.2 • 10" 4 0.05 1.85 9.6 4096 3 
460 0.189 3.6 • 10~ 3 2.5 • KT 5 3.07 • 10~ 3 1.45 • 10~ 3 0.083 1.85 9.9 2048 3 
320 0.187 3.5 ■ KT 3 0.5 ■ KT 4 6.14 ■ 10" 3 2.45 • 10" 3 0.12 1.85 10 1024 3 



Table 1. Parameters of the numerical simulations. 5Ra = y/^u ims L/ v: Taylor- 
Reynolds number, u rms : root-mean-square velocity, e^: mean kinetic energy dissipation 
rate, v. kinematic viscosity, dx: grid-spacing, rj — (y 3 /e^) 1 / 4 : Kolmogorov dissipation 
length scale, t v = (ly/e^) 1 ^ 2 : Kolmogorov time scale, L = (2/3E) 3 ^ 2 /e] i : integral scale, 
Tl = L/u Ims : large-eddy turnover time, iV 3 : number of collocation points. 

two in the case of 3?a — 730). Parameters of these high- Reynolds number simulations 
are given in Table [TJ 

Fig. [2] shows the application of the rescaling formula (IE]) to the 2nd, 4th, 6th and 
10th order structure functions obtained from a Navier-Stokes simulation with 2048 3 grid 
points and parameters as described in Table (II]). We choose this data set because it 
contains ten large-eddy turn-over times and thus provides reliable statistics for high- 
order structure functions. In each sub-figure of Fig. EJ the unsealed structure functions 
are shown in the lower part (solid lines) and the structure functions rescaled according 
to eqn.(|H]) are show on top (dashed lines). Note that the original structure functions 
are shifted for clarity. The effect of the rescaling-transformation (jHJ) is two-fold: first, 
the amplitudes from the dissipative scales up to the integral scales are now very similar. 
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Figure 1. Comparison of structure functions: longitudinal (black), transverse (blue), 
rescaled longitudinal (red) and the one using the differential relation © (green). 
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Figure 2. Structure functions of different order (5R> = 460). dashed lines: rescaled 
abscissa according to (6); solid lines: original (shifted vertically for clarity), Inset: 
logarithmic derivatives, black diamonds: rescaled transverse structure functions, 
vertical lines indicate the flat region of the third order structure function. 



In addition, the range of scales over which a power-law behavior is observed is now 
identical for the two SFs, although the scaling exponents differ slightly. This is evidence 
in the inset of each sub-figure in Fig. |2] where the logarithmic derivatives of the structure 
functions with respect to scale have been plotted, and the vertical lines mark the scaling 
interval. Note, that both effects could not be achieved by an order-independent fixed 
rescaling factor of 3/2. 

For this data set the transverse increments are more intermittent than the rescaled 
ones. In order to address the question whether this difference in scaling might depend on 
Reynolds number, we show in Fig. |3]the logarithmic derivative of the 8th order structure 
function for three different simulations. With increasing Reynolds number the inertial 
range increases but the scaling exponent (value of the plateau) remains the same. We 
measure a value of approx. 2 for the transverse functions and 2.21 for the longitudinal 
ones. For comparison we included data from randomized snapshots originating from the 
simulation with 3?a = 730. In detail, for each Fourier mode we change randomly the 
phase while preserving its amplitude and incompressibility of the flow. This preserves 
the energy spectrum but destroys the structure of the flow (energy cascade, coherent 
structures ...). The randomized longitudinal and transverse structure functions exhibit 
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Figure 3. Logarithmic derivative of longitudinal and transverse structure functions 
for three different Reynolds numbers and in the case of a randomized velocity field. 

the same scaling exponent, now close to the trivial 8/3 value. 

A general assumption is that possible remaining large scale anisotropy in the small 
scales is expected the decrease with Reynolds number. We remark that Biferale, Lanotte 
and Toschi p3] showed that the differences in the high-order exponents remain even if 
measured in the purely isotropic sector. That the curves in Fig. [3] fall on top of each 
other within the inertial range of scales is an indication that the observed differences 
in the scaling exponents are not due to large scale anisotropies. This indicates that the 
former observed differences have to be attributed to the specific small scale structures 
of the flow. 



3. Implications for longitudinal and transverse PDFs 

Since the rescaling property has the effect to make the longitudinal and transverse 
structure functions fall nearly on top of each other, we want to understand the effect 
of the rescaling transformation on the probability density functions (PDFs). In this 
subsection we try to map longitudinal PDFs to transverse ones using the rescaling 
property expressed through eqn. (JHJ). The rescaling transformations were derived for 
even order structure functions. Thus in the following, we disregard skewness effects and 
consider only the symmetric part of the PDFs. We first approximate the numerically 
obtained longitudinal PDFs with a log-normal distribution using the expression given 
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in Yakhot 

r 'In. 



P L (Au,r) = — -== / e x exp 

TiAuvinr 



oo 

2 



461nr 



cfe (10) 



for which a fit is obtained with the values a = 0.383 and b = 0.0166 [22]. In Fig. H]the 
numerically obtained PDF and the fit Pi(Au,r) are shown for several spatial scales. 
We apply the inverse Mellin transform 



I rioo 

P L (Au,r) = — dnS(n,r)(Au) 
Am J_ ioo 



with S(n,r) = A(ri)r^ n \ and we follow the procedure in [22] which fixes the amplitude 
by going to the Gaussian limit for large spatial differences: 

on/2 poo 

A(n) = (n-l)\\ = —= / e~ x2 x n dx 



Now a mapping from the longitudinal PDFs to the transverse PDFs is obtained by 
inserting the rescaling relation (jSJ) in the expression for the structure functions: 

i pioo 

P T (Au,r) = — dnA{n){C{n)rf {n \Au)- n (11) 
Au J_ ioo 

where C{n) = yj^l+ij as in ec l n - ©• 

Since the ansatz S(n,r) = A(n)r^ n ^ does not contain a cutoff at integral and at 
dissipation scales, this cutoff is inserted in the rescaling function C(n). In both regions 
outside the inertial range smooth behavior is expected. On scales close to the integral 
range, Gaussian behavior for both longitudinal and transverse increments is expected 
and thus no rescaling is necessary. This justifies to choose C(n) to be constant for 
n < 0. The cutoff at the dissipation scale is achieved by choosing C{n) to be constant 
for n > 6. The precise value of the chosen n is dependent on the actual Reynolds number 
and the effect of choosing a different bound allows for a widening of the transformed 
PDF. This reflects the fact the Reynolds number has a similar effect on the width of 
the PDF. Evaluating the integral ( ITT]) using a saddle point approximation, we obtain 
a mapping from the log-normal fit of the longitudinal PDF Pi(Au,r) to a new PDF 
Pt{Au, r) which is compared to the numerically obtained data in Fig. H]for increments 
ranging from the near dissipation range to integrals scales. One may observe that the 
agreement is especially remarkable in the inertial range (r = IO677 and r = 212^.) 
We do not expect perfect agreement for all scales since in that case there would be 
no room for differences between longitudinal and transverse structure functions. Thus 
the discrepancy for r = 21r] and r = 427/ just represents the missing contributions of 
the pressure term. Therefore, this method of mapping the PDFs is also a promising 
candidate for applications like PDF modeling of turbulent flows (see Pope [23] and 
references therein). 
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Figure 4. Longitudinal PDF (lower points) and fit with log-normal distribution (green 
line); transverse PDF (upper points) and mapped distribution (red line). Shown are 
the PDFs for different increments ranging from the near dissipation range to integrals 
scales. (3?a = 460, w denotes the dissipation and L the integral scale.) 

4. Conclusions and outlook 

In this paper, we have suggested a new way of analyzing experimental and numerical 
data for longitudinal and transverse structure functions in Eulerian data of a turbulent 
velocity field. This procedure yields a mapping between the longitudinal and transverse 
scales, which provides consistent reference point in the identification of the inertial 
ranges of scales of turbulent flows. In addition, the derived scale correspondence allows 
for a direct mapping of the full probability density of transverse and longitudinal struc- 
ture functions. This may be of much practical interest as the distributions carry a more 
complete information that than a subset of their moments. The gap of longitudinal 
and transverse structure function exponents seems not to depend on Reynolds number 
but on small scale structure of the flow. The proposed mapping may help clarify the 
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role played by the pressure terms. Future work will be devoted to the analysis of other 
turbulent systems like magnetohydrodynamics, where the addition of the magnetic pres- 
sure term poses an interesting comparison. 
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